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Organization of strongly interacting directed polymer liquids 
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The impact of impenetrable obstacles on the energetics and equilibrium structure of strongly 
repulsive directed polymers is investigated. As a result of the strong interactions, regions of severe 
polymer depletion and excess are found in the vicinity of the obstacle, and the associated free- 
energy cost is found to scale quadratically with the average polymer density. The polymer-polymer 
interactions are accounted for via a sequence of transformations: from the 3D line liquid to a 2D 
fluid of Bose particles to a 2D composite fermion fluid and, finally, to a 2D one-component plasma. 
The results presented here are applicable to a range of systems consisting of noncrossing directed 
lines. 
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Many investigations of (d+l)-dimensional systems of 
directed lines hinge on an analogy between the statistical 
mechanics of these systems and the quantum mechan- 
ics of d-dimensional systems of nonrelativistic bosonic 
particles pQ. This analogy has proven useful in work 
on directed line liquids, such as those consisting of 
directed polymers (either under tension [5] or within a 
nematic solvent [3J) as well as fluctuating line defects in 
symmetry-broken phases, such as vortex lines in type-II 
superconductors [4]. 

In this Letter, we focus on the specific example of 
a liquid of many thin thermally fluctuating polymers 
under tension, which interact with one another via a 
strong short-range repulsion. We represent this liquid in 
terms of (2+l)-dimensional directed lines, and we treat 
the polymer repulsion as a restriction on the possible 
configurations, prohibiting the intersection of any two 
lines. Central to the present work is the non-perturbative 
treatment of this nonintersection restriction, which we 
accomplish via the well-known Chern-Simons transmu- 
tation of quantum statistics, from Bose to Fermi. Using 
quantum many-body theory, we address two fundamental 
questions regarding this system, viz., "What are the free- 
energy cost and equilibrium polymer density induced 
by an impenetrable inclusion, i.e., a region in which 
the polymer density is forced to be zero?" We note 
that the free-energy cost of the inclusion corresponds 
to the probability that the unconstrained system would 
spontaneously satisfy the constraint. For simplicity, we 
study thin, flat inclusions within a planar slice perpen- 
dicular to the polymer direction z, and we determine 
the equilibrium polymer density within that slice. More 
generally, we determine the free energy and equilibrium 
density in the presence of spatially extended constraints, 
such as a ring, positioned in the slice, threaded by a fixed 
number of polymers, as shown in Fig. [I] The inclusion is 
precisely the case for which no polymers thread the ring. 
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FIG. 1. [Color online) Sequence of transformations of a 
liquid of directed polymers subject to a spatially extended 
constraint, shown as a ring (thick line), threaded by a fixed 
number of polymers. The sequence is used to obtain the 
free energy and polymer density of the constrained system. 
Ultimately, the directed polymer liquid is transformed to a 
two-dimensional one-component plasma, whose the Coulomb 
repulsion results from integrating out the long noncrossing 
chains. 



We emphasize that the dominant free-energy cost of 
such constraints is the result of the many-body effect of 
the crowding of the polymers, an effect absent for free 
directed lines. 

The transformation between directed lines and bosons 
enables one to harness powerful techniques from quan- 
tum many-body physics for application to directed line 
liquids. The first such application was due to de Gennes, 
who invoked it to calculate the structure factor for 
a (l+l)-dimensional hard-core directed line liquid [2J. 
A recent elaboration on these results by Rocklin et 
al. yielded the free energy and equilibrium polymer den- 
sity in the presence of pins and other impenetrable con- 
straints for the same system [5]. In essence, de Gennes's 
strategy separates into two major steps. First, the 
classical (l+l)-dimensional noncrossing directed lines are 
mapped into quantum hard-core point bosons moving in 
one dimension pQ. Second, these hard-core bosons are 
mapped, exactly, to free fermions [6]. The virtue of this 
sequence of transformations, which we are here seeking 
to emulate for (2+l)-dimensional systems of noncrossing 
directed lines, is that it reduces the strongly interacting 
directed line liquid to an (exactly solvable) system of 
free fermions. We accomplish a similar reduction by 
mapping our (2+l)-dimensional directed line liquid to a 
two-dimensional fluid of Bose particles, and then make 
the Chern-Simons transmutation from Bose to Fermi 
statistics [7J. 

We address a model of the line liquid that consists of 
N paths that are directed along the z-axis from z = to 
z = L via a line tension r, which penalizes deflections of 
the paths. We assume no path switchbacks, in which case 
the (Cartesian) (x, y)-coordinates of the N paths R(z) = 
{r n (z)}n = i are single- valued functions of z. We study 
interactions that are short-ranged and so strong that they 
prohibit configurations in which any pair (n, n') of lines 
intersect, i.e., in which r„(z) = r n >(z) for any value of 
z. The statistical weight for noncrossing configurations 
is then proportional to exp {—Q/T), where 
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and the temperature T is measured in units of the Boltz- 
mann constant. Correspondingly, the partition function 
Z is given by the functional integral Z = J[-DR(-)]e~ e / T , 
taken over all noncrossing paths [8]. 

To implement the mapping to quantum many-body 
physics we note that Z = ($ F \e' 0Ln \^ 1 )} see, e.g., 
Ref. [T|. The quantum Hamiltonian, 



H 



1 



A' 



n— 1 



(2) 



describes two-dimensional Bose particles of mass m under 
the correspondence (r, T) <H- (m,K). In addition, the 



quantum amplitudes (R^ 1 ^) reflect the initial and final 
a priori classical end-point distributions. The quantum 
states |R) are symmetrized products of N single-particle 
position eigenstates |r), appropriate for indistinguishable 
polymers. Observe that the limit L — > oo is also the low- 
temperature limit /3 — > oo, for which the quantum ground 
state dominates Z. 

Our focus is on the increase of the free-energy AT 
(= T c — T u ) from T u (its value in the unconstrained 
directed line liquid) to T c (its value in the liquid with 
a spatially extended constraint). We compute AT using 
the ratio between the unconstrained partition function Z 
and the partition function Z c , which includes the effects 
of constraint. If, as in the case we consider, the constraint 
requires that exactly Q directed lines pass through some 
two-dimensional region T> perpendicular to the director z 
and far from the ends z — 0, L, then, in the limit L — > oo, 
AT may be computed using ground-state dominance as: 



AT = -Tin— = -Tin 
Z 



dR |* fc (R)| 5 



(3) 



Here, C indicates the constraint on integration that 
exactly Q of the N coordinates {r„} lie within T>. The 
ground-state wavefunction in Eq. ^ *&(R) (= (R|GS)) 
is the lowest-energy (particle-exchange symmetric) solu- 
tion of the energy eigenproblem H^b = ESSft,- Thus, we 
are confronted with the task of finding an expression for 

l* h | 2 - 

The restrictions r„(z) ^ r n >(z) on the path integrals 
for the line-liquid partition functions Z and Z c demand 
a nonperturbative treatment. Following the essence of 
de Gennes's [2J strategy, we transmute the statistics of 
the quantum fluid from Bose to Fermi [7], thus arriving 
at a description in terms of a many-fermion wavefunction 
^ f that necessarily has nodes for any two coincident 
particles, i.e., \P/(ri, . . . ,r, . . . ,r, . . . ,rjy) = 0. This 
transmutation is accomplished via the well-known [9 
singular gauge transformation of Chern-Simons theory: 

* b -)• */ = *&exp(i<^„'<n0™,«'): where n , n > = 
tan" 1 [(y„ — y n >)/(x n — x n >)] is the polar angle between 
particles n and n'. Then, under the exchange of particles 
n and n', n>n > — > 9 n ,ri ± it, and, provided <f> is an odd 
integer, $>f is antisymmetric, i.e., fermionic. As ^b is 
the ground state of H, Eq. (pi, ^ f obeys the transmuted 
energy eigenproblem [9] T-L'^f = E'^f, where 
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W has the same energy spectrum as W, so for any 
eigenstate, including the ground state, E' = E. 



If, in Eq. (4a), A„(r n ) were to depend only on r„ and 



not on any {r n /(^ n )}, the Hamiltonian 7i' would describe 
independent particles in a magnetic field V x A. However, 
as shown in Eq. (4b), A„ does depend on {r n /(^„)}, 
which implies (nonlocal) interactions between all parti- 
cles. Equation ( [4b] ) states that {A„} describes particles 
that have <p quanta of fictitious magnetic flux attached 
(i.e., localized at the position of every particle, each 
quantum carrying h/q flux, where q is a fictitious charge) . 
Thus, Ti! describes composite fermions — composed of 
particles obeying Fermi statistics and flux tubes (TU1 HJ • 
This transmutation of statistics is the Hamiltonian form 
of Chern-Simons theory [12] . 

Although this formulation effectively incorporates the 
hard-core restriction, it remains intractable, and thus we 
are led to take advantage of a natural approximation (see, 
e.g., Refs. [IQlUT]), which we now describe. Instead of 
having an odd number of flux tubes attached to each 
particle, we smear the magnetic field associated with one 
flux tube per particle uniformly over the area of the 
system, and gauge-transform away the remaining flux 
tubes. In this so-called Average Field Approximation 
(AFA), the fermions are non-interacting and subject to a 
homogeneous magnetic field VxA = Bz, corresponding 
to one quantum of magnetic flux per particle. [In the 
symmetric gauge, A = B(y, —x, 0)/2.] In the language 
of directed lines, the magnetic field is equal to the 
number of lines per unit area po times a quantum of 
flux (i.e., B — p h/q). In this magnetic field, the many- 
body ground-state has energy E = NBhq/2m and 
the wave-function Slater determinant takes the Vander- 
monde form |13) 



#?(R) ex e 



-££Ul^l 2 /« 2 



n ( w " 



(5) 
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where w n = x n + iy n . Within this approximation, the 
quantum ground-state energy and wave-function describe 
a long, noncrossing, directed line liquid. In particular, for 
its free energy per unit volume g (= (Q)/AL), we obtain 
the result 



g = TiT 2 pl/T. 



(6) 



As Eq. (H)]) gives an expression for the ground-state 
wave-function, we proceed with our goal of calculating 
the increase in the free energy AT by using the integral in 
Eq. ([3]). We thus make the approximation \^b\ 2 ~ [*/ 1 2 > 
and compute AT = -Tin J c dXexp (-[/(R)/T), where 



U(R 



N 



T ln|^| 2 = 2]Tln|r n - 

l<n<ro'<JV 



r n ,\-7rp J2\rn\ 2 . (7) 



Interpreted as a potential energy, U describes a two- 
dimensional one-component plasma (2DOCP) [T51 [T5] . 
In the general case of the plasma, particles of (e.g., 




FIG. 2. (Color online) Two-dimensional one-component 
plasma (2DOCP) subject to a spatially extended constraint. 
The electrostatic approximation predicts a gap in the mo- 
bile charge density and an accumulated surface charge for 
cases (A), for which the excess mobile charge resides outside 
the ring, and (B), for which the excess mobile charge is 
confined inside the ring. This approximation also predicts 
discontinuities in the mobile-charge density, which in the 
exact solution are smeared out due to thermal fluctuations 
of the mobile charges around the minimum energy state. The 
exact solution leads to the density profiles p(r)/po shown 
for (Q - Qo)/V2Qo = ±4 (thick) and (thin), where -Q 
and Qq (= -ira 2 po) are the mobile and background charges 
inside the ring, respectively [141 115] . The inset shows the 
rescaled energy cost (4/Qq) AT [top: case (A), bottom: case 
(B)] [TB]. The curves are valid for macroscopic values of 
the charge deficiency, i.e., for not-too-small values of the 
argument, when AT is dominated by electrostatic energy. In 
the text, we establish that these profiles corresponds to the 
in-plane distribution of noncrossing directed lines due to a 
planar spatially extended constraint. 



negative) charge — e inhabit a uniform background that 
maintains overall charge neutrality, and they interact via 
the two-dimensional Coulomb repulsion: e 2 ln|r„ — r„<|. 
This U/T corresponds to a specific value, viz. 2, of the 
plasma coupling constant e 2 /T. Applied to the plasma, 
the measure J„ dH in AT requires exactly Q mobile 
charges to occupy the region of constraint T>. 

This completes our reduction of the three-dimensional 
classical liquid of noncrossing directed lines to the clas- 
sical two-dimensional one-component plasma (2DOCP); 
this reduction enables us to compute physical properties 
such as the free-energy cost and the equilibrium density 
profiles associated with the imposition of spatially ex- 
tended constraints. The long-range planar interactions 
within the 2DOCP arise from the short-range interac- 
tions between the fluctuating directed lines, integrated 
over their length. To calculate the free energy and the 



equilibrium density of the directed line liquid in the 
presence of a ring constraint, we analyze the 2D0CP, 
first at the level of classical electrostatics, and then 
allowing for fluctuations. Note that the constraint results 
in regions of excess mobile charge and regions partially 
depleted of mobile charge. To minimize the electrostatic 
energy, any excess mobile charge is forced up against 
the boundary of T>. Moreover, on the other side of the 
boundary, a region fully depleted of mobile charge (i.e., 
a gap) opens up, out to a radius within which the net 
charge is zero [17] , 

Similarly, at the level of electrostatic energy mini- 
mization one can readily determine the energy of the 
charge distribution and, hence, the free-energy cost of the 
constraint AJ r (Q, Q ), in terms of the mobile charge — Q 
and the background charge within T>, i.e., Qq (= ira 2 po), 
both in units of e [16] ; see Fig. [2] For the special case 
Q = (i.e., T> empty, corresponding to an inclusion) 
the free-energy cost has a simple form: A J 7 = Qg/4 = 
7r 2 a 4 pQ/4. This electrostatic result is significant for the 
noncrossing-directed-line liquid, as it differs qualitatively 
from the case of a noninteracting directed line liquid, for 
which a particle inclusion incurs a free-energy cost linear 
in a 2 p a . 

To improve upon the electrostatic approximation, we 
take into account the effect of thermal fluctuations on the 
polymer density profile. We rely on the exact solution 
of the 2DOCP with the appropriate plasma coupling 
constant, i.e., 2; see Refs. Q3HH3 UM ES2 • In the limit 
a ^ 1/a/PO; the exact density profile outside the region 
of constraint depends only on a, p and Q through the 
combination (Q — Qo)/V^Qo [II]- The layer of excess 
mobile charge on one side forms an electrical double layer 
of thickness of order l/^/p^; see Fig.pt and Ref. |14) . The 
region partially depleted of mobile charge does develop 
a soft gap, in which the charge is small but nonzero. 
The mobile charge density profile progresses smoothly, 
according to a qualitatively error-function-like curve, 
through the boundary region, rapidly approaching the 
value that exactly compensates the background charge 
density po; see Fig. [2b and Ref. [T4"] . The mobile charge 
density profile for the depleted side of the constraint 
applies to both cases, Q > Qo and Q < Qq, and similarly 
for the excess-charge side. For Q small relative to Qo, 
the remaining mobile charge in T> forms a droplet whose 
shape is essentially the density profile of a system of 
electrons that fill the lowest Landau level: a flat central 
profile and a decay into the soft gap; see, e.g., Ref. [2T] , 
Thus, we have established that when some fixed portion 
of the lines of a directed line liquid is constrained to 
thread T>, the equilibrium density profile in the slice 
containing T> is that of the correspondingly constrained 
2DOCP. 

Recapping our strategy, we progressed from a three- 
dimensional liquid of thermally fluctuating lines, to a 
two-dimensional quantum many-boson fluid, to a two- 



dimensional quantum many-fermion fluid coupled to a 
Chern-Simons gauge field, which we treated in the Av- 
erage Field Approximation to obtain the filled lowest 
Landau level picture. The phenomenology of a lowest 
Landau level filled with noninteracting fermions is well 
studied, and suggests various analogous phenomena for 
the corresponding hard-core boson fluid. However, as the 
AFA is an approximation, these analogous phenomena 
may be artifacts of the approximation, and we now use 
physical intuition to identify any such artifacts. For 
example, the quantum Hall effect suggested by the AFA 
is one such artifact: the boson fluid does not have 
broken time-reversal symmetry, and therefore shows no 
Hall effect [25H27] , A second artifact is suggested by 
the incompressibility of the filled lowest Landau level, 
which would incorrectly imply the incompressibility of 
the boson fluid. However, by reinstating the inter- 
particle interactions A(r) — A a (r), e.g., via the Random 
Phase Approximation, the compressibility of the boson 
fluid is restored, as shown in Refs. [55] [25]. In the context 
of the directed line liquid, the thermodynamic (areal) 
compressibility k is defined in terms of the free energy 
density g via n^ 1 = p\ d 2 gjdp\, and thus, using Eq. (tfjj), 
we obtain k = t/2ttT 2 Pq. A particularly noteworthy 
consequence of the residual interactions A(r) — A a (r) is 
their ability to renormalize the effective plasma coupling 
constant e 2 /T in the plasma analogy away from the 
exactly solvable case, viz., 2. Nevertheless, we expect 
the general picture presented here, of the energetics and 
structure of the directed line liquid in the presence of 
spatially extended constraints, to hold. 
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